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We present some results pertaining quantities which are regarded as 
good indicators of the pseudo-critical temperatures for the deconfinement 
and partial chiral restoration transitions using a polynomial form for the 
Polyakov potential part and an extended version of the NJL model which 
includes 6 and 8 quark interaction terms. Some comparisons with results 
from the lattice formulation of QCD (IQCD) are performed and results for 
the location of the critical endoint in the phase diagram are also presented. 
It is shown that the comparison with IQCD results favors a moderately 
strong OZI-violating 8 quark interaction term. 

PACS numbers: ll.lO.Wx; ll.30.Rd; 11.30.Qc 



1. Introduction 

Low energy effective models are a valuable tool for the study of the 
low-energy thermodynamic properties of strongly interacting matter. Their 
importance is even greater in the case for finite chemical potential where 
the sign problem afflicts the more fundamental approach of using IQCD. 

Two of the main features of the non-perturbative regime of QCD are 
the spontaneous breaking of chiral symmetry and confinement. For both a 
transition is expected to occur at finite temperature and/or chemical po- 
tential: partial restoration of chiral symmetry and deconfinement (with the 
spontaneous breaking of the Z{Nc) centre symmetry). 

The NJL model [T] incorporates by construction a mechanism for the 
dynamical breaking of chiral symmetry and shares with QCD its global 
symmetries. The unwanted axial symmetry can be removed by the inclusion 
of a 't Hooft determinantal [21 [3l |3] term but this however introduces a 
vacuum stability problem when considering the light quark sector (u, d and 
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s quarks). This issue can be solved by further extending the model to 
include eight-quark interaction terms O |6] . 

A partial inclusion of the gluonic degrees of freedom can be done through 
the consideration of a static background homogeneous temporal gauge glu- 
onic field coupled to the quark fields through the covariant derivative |7|. 
The Polyakov loop can then be used as an approximate order parame- 
ter for deconfinement (it vanishes for the confined phase). An additional 
temperature-dependent pure "gluonic" term, the Polyakov potential, U, for 
which several forms have been proposed, must also be introduced to drive 
the transition. Due to space constraints here we will only present results 
referring to the polynomial form as proposed in [8] (furthermore the results 
are qualitatively similar with the other forms tested). 

The specification of an non-renormalizable effective model is not com- 
plete without the choice of a regularization procedure. In the present study 
we will use a Pauli-Villars regulator with two subtractions in the integrand 
[S] . In |1U] we have shown that using this we can obtain the expected high- 
temperature asymptotic behaviour for the PNJL model while consistently 
using the same regularization procedure for both the vacuum and medium 
contributions for the relevant integrals . 

The thermodynamic potential is given by (for more details see \10\ 111)): 

n{Mf,T,fi,cP,^) = 
=^ (^4Gh} + KKhdhs + ^ {h})^ + ^g2h) 

+^ E {J^l{MlT,^^,<i,A) + C{T,^x))+U{ct),'4>,T). (1) 

f=u,d,s 

Here G, k, gi, §2 refer to the coupling strengths of the 4q (NJL term), 6q 
('t Hooft determinant) and OZI-violating and non-violating 8q interactions 
respectively. Mj and hf {f = u,d, s) refer to the dynamical masses and the 
condensates which must satisfy the conditions (m j refers to current masses) : 

m„ - M„ = Ghu + %hdhs + ^Khj + ^hl 
md-Md = Ghd + %Khs + S^hdhj + f /i3 (2) 
ms-Ms = Ghs + §Khd + ^hshj + f /if. 

J_i and C (a mass independent, but T, // dependent, term which must be 
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added to ensure thermodynamic consistency |10l [TT] ) are given by: 

/>oo 

rr (m2) = - / dp pa8p^ {Em - Eo) 
Jo 

J^f {M',T, p, <A, 0) = r dp 8p'pA (2 {Em -Eo) + Tln 

Jo \ riqonqo 

c(r.M) = -|r*/f==^^^V (3) 



V P 

Here E is the energy and Uq/n-q refer to the occupation numbers (the 
subscript refers to the occupation number for massless particles). The 
extension to include the Polyakov loop in the J_i integrals can be done 
straightforwardly by noticing that its phase enters the action as an imagi- 
nary chemical potential and modifying the occupation numbers accordingly. 
The regulator was chosen to be pA = 1 — (l — Iv^Sm'^) exp {K^Om'^^ 

2. Transition signals 

Several quantities are expected to undergo a rapid change during the 
transitions acting therefore clS cl SI gnal for the determination of the pseudo 
critical temperatures (note that with physical quark current masses both the 
chiral and the Z{Nc) centre symmetries are only approximate; they are only 
exact in the chiral and quenched limit respectively) . Besides a rapid change 
in the condensates and Polyakov loop one also expects a distinct signal 
coming from the chiral, quark number and Polyakov loop susceptibilities: 



1 / 5217 



r2 V 92 



^=0/ ^"""^ ~ dpi ' 



- 1 2^-^ (4) 

The use of the dual chiral condensates has been suggested as a possi- 
ble additional indicator for deconfinement |12J. These correspond to the 
Fourier transform with respect to an angle which is introduced through the 
consideration of twisted temporal boundary conditions: 

= J ^ -^h,{a)da, (5) 

where hi{a) is given by ^ = with p ^ p + iTa, (with — vr < a < vr). 
Results for the temperature dependence of the dressed Polyakov loop (which 
corresponds to setting j = 1 in the above) will also be presented here. 
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3. Results 

As we reported in [6] it is possible to obtain the same low energy scalar 
and pseudoscalar mesonic spectra throughout a large range of values for 
the OZI-violating 8q coupling strength by simultaneously lowering G when 
increasing gi (except for the a meson mass which decreases for higher gi). 
Here we will focus therefore in the impact of this choice in the transitions 
behaviour. Another parameter, Tq, coming from the polynomial form of 
U has also to be taken into account (it sets the temperature scale for the 
temperature dependence of U). 

As can be seen in Fig. 1 the chiral condensates, the Polyakov loop and 
the dual chiral condensates undergo a rapid change in a narrow temperature 
range signaling the transitions. The change is more rapid and occurs at a 
lower temperature for stronger 8q interactions. The behavior of the dressed 
Polyakov loop closely follows that of the chiral condensate and the zeros of 
the second derivative with respect to the temperature are in fact very close 
to each other (a few MeV apart at most). The temperature dependence of 
the susceptibilities is also strongly influenced by the value of gi as can be 
seen in Fig. 2: the peaks get sharper and higher with increasing gi 




(a) (b) 

Fig. 1. Temperature dependence of: a) the Polyakov loop (full lines starting at zero) 
and chiral condensates divided by the peak value (light in full lines and strange 
with dashed lines); b) the dressed Polyakov loop (light in full lines and strange in 
dashed lines). The different sets correspond to different values for gi in the PNJL 
model, 1^\to=i9o MoV 5i ^ 10"'^ = 1, 5, 6.5, 8 GeV^^ (for the chiral condensates and 
dressed Polyakov loop lower lines correspond to higher gi whereas for the Polyakov 
loop it is the opposite). 

A comparison with IQCD resuhs (see Fig. 3 a) and b)) for several quan- 
tities such as the interaction measure and the subtracted chiral condensate 
gives a reasonable result. The deviation from the IQCD result in what refers 
to the subtracted chiral condensate is probably related to the slow chiral 
restoration of the strange quark (Fig. la). The position of the critical end- 
point of QCD in the phase diagram is also very affected by the value of gi 
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(d) (e) 

Fig. 2. Temperature dependence of susceptibilities in the PNJL model iJA with Tq = 
190 MeV) for different values of gi at vanishing chemical potential. Top row from 
left to right: light quark chiral susceptibility, strange quark chiral susceptibility 
and Polyakov loop susceptibility. Bottom row: temperature derivative of the light 
and strange quark number susceptibilities. In all of these panels the sharper peaks 
correspond to stronger OZI- violating 8q interactions (same sets as in Fig. 1). 



as can be seen in Fig. 3c. The Tq parameter from U. on the other hand only 
seems to affect the temperature at which it occurs and not the chemical 
potential. 
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(a) (b) (c) 

Fig. 3. Interaction measure {gi = 6000 GeV"^ Tq = 0.190, 0.210, 0.230 GeV; 
higher Tq results in a lower and more shifted to the right peak) and the subtracted 
chiral condensate {Tq — 210GeV for comparison we also show in dashed lines the 
result obtained with the logarithmic form oiU from [14]) compared to IQCD results 
(data taken from jl3l ). Temperature (positive slope lines) and chemical potential 
(negative slope lines) of the CEP as a function of gi for Tq — 190, 210, 250 MeV 
(higher Tq corresponds to thicker lines) as well as for the NJL model. 
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